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ABSTRACT 



We present an extensive search for a general class of nipped SU(5) models built within 
the free fermionic formulation of the heterotic string. We describe a set of algorithms which 
constitute the basis for a computer program capable of generating systematically the mass- 
less spectrum and the superpotential of all possible models within the class we consider. 
Our search through the huge parameter space to be explored is simplified considerably 
by the constraint of N = 1 spacetime super symmetry and the need for extra Q, Q rep- 
resentations beyond the standard ones in order to possibly achieve string gauge coupling 
unification at scales of 0(1O 18 GeV). Our results are remarkably simple and evidence the 
large degree of redundancy in this kind of constructions. We find one model with gauge 
group SU(5) x U(l)~ x SO(10)h x SU{A)h x £7(1) 5 and fairly acceptable phenomenolog- 
ical properties. We study the D- and F-flatness constraints and the symmetry breaking 
pattern in this model and conclude that string gauge coupling unification is quite possible. 
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1. Introduction 

String model-building has provided a new and profuse source of ideas to overcome 
many of the weaknesses of traditional unified models. Indeed, any string-derived model 
possesses a definite gauge group whose gauge couplings unify at a calculable high-energy 
scale, and its matter representations and their gauge and Yukawa interactions are com- 
pletely determined. In practice this means that one faces a rather constrained problem 
with little or no room for ways out of potentially phenomenologically disastrous situations. 
Fortunately (or unfortunately) there is a very large number of models to choose from and 
one then hopes that an educated path through the morass of possibilities might lead to 
a model which describes the features of the low-energy world and also predicts new ob- 
servable phenomena. The purpose of this paper is to pursue one such path through a 
(hopefully sizeable) portion of this space. 

There are several classes of constraints that one can apply to the universe of possibil- 
ities to reduce the sample that needs to be considered. The most basic ones are: N = 1 
spacetime supersymmetry and a gauge group below the string scale that either includes 
the standard model gauge group or embeds it but can be dynamically broken down to it at 
a lower energy scale. There is also an array of phenomenological constraints that need to 
be satisfied by the successful model: the low-energy spectrum must contain the three gen- 
erations of quarks and leptons, and the very accurately measured values of sin 2 Q Wl a e , and 
«3 must be reproduced. The fulfillment of both these constraints depends on the details of 
the path from the Planck scale down to low energies, and are therefore harder to enforce 
given our present lack of understanding of the mechanics of supersymmetry breaking. At 
a deeper level of detail we must reproduce the observed values of the quark and lepton 
masses. Any string-derived model which satisfies all these constraints is a candidate for the 
correct theory of all particles and interactions, and can be tested by e.g., their predictions 
for the top-quark, neutrino, higgs, and sparticle masses. 

From the string theory point of view there are some organizing principles. It is widely 
believed that all string models can be thought of as points in a parameter space called the 
moduli space. Each particular model having fixed values of these parameters, although 
classes of models with continously connected values of the parameters are very common. 
The hope is that eventually a string theory principle will be found which will somehow 
select the energetically most favorable point in this space. In practice, this space has 
been parceled up into several different 'formulations' of string theory [pj which describe 
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(not necessarily non-overlapping) subsets of models. For model-building purposes we have 
chosen to explore (a portion of) this parameter space with the visor of the so-called free 
fermionic formulation of the heterotic string |^,[^[|,|5|,[)],|7| . This formulation can be described 
as a set of mo del- building rules which is amenable to systematic algebraic and symbolic 
manipulation and therefore can be translated into a computer code. The description 
and implementation of such code is one of the objectives of this paper. A further bonus of 
working within this formulation is the ease with which the terms in the low-energy effective 
action can be calculated ||, such as the cubic and higher-order superpotential couplings. 

String theory also provides a correlation between the gauge group and its allowed 
matter representations. First of all, these always come in anomaly-free sets (except for a 
possible pseudo-anomalous Ua(X) subgroup). The four-dimensional gauge group is a reflec- 
tion of the algebra of two-dimensional world-sheet currents called the Kac-Moody algebra 
||. These algebras can be realized at integer values of the 'level' k = 1, 2, . . ., all of which 
represent the same gauge group. The allowed matter representations depend on the chosen 
level which basically determines a cap on the dimensionality of the allowed representations 
at that level [pT0|JTT|] . By far the simplest and most common realizations are level one. 
These allow a very limited set of matter representations which do not include the adjoint 
representation. This result eliminates from further consideration any of the traditional 
unified groups (e.g., SU(5), SO (10), Eq). Higher-level realizations are possible (and they 
allow the adjoint representation) although in practice these are harder to construct |IT] . 
Besides, the unified groups that would become phenomenologically viable require high lev- 
els of the Kac-Moody algebra to accommodate the large representations which appear in 



traditional model-building with these gauge groups [lQ|. Limiting ourselves to level-one 
Kac-Moody algebras we have two classes of models to consider: the class of 'flipped'-like 
SU(n) models (of which flipped SU(5) [p~3| , pT^JT5|] is the simplest) and related models pf , 



or models which contain explicitly the standard model gauge group [|17l , |18| , |T9[1 and there- 
fore require no further symmetry breaking. In this paper we explore the simplest unified 
models which can be constructed using the simplest Kac-Moody algebras, i.e., the class of 
flipped SU (5) string models. 

There already exists one string-derived realization of a flipped SU(5) model in this 
formulation, the so-called 'revamped' flipped SU(5) model p0| , pl| , |22[] . This model has 



been explored in detail in the literature P^,p5|,p6|,[27|,P8|, 29| and several of its interesting 



1 Flipped SU(5) models have also been constructed in other formulations [23]. 
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features have been highlighted, such as the natural apperance of a hierarchical fermion mass 
spectrum P^ , p7[ , adequate low-energy higgs sector pip7|] , acceptable proton decay rate 



|26|1 , the existence of bound states of fractionally charged hidden sector particles called 
cryptons |25[], etc. Despite all these attractive features, gauge coupling unification at a 
scale of O(10 18 GeV) as predicted in this model [ |3"U| , |3"T| ] , is probably not possible 



or 



conversely, the low-energy values of sin 2 9 W and as are probably not reproducible. In any 
event, one must consider these statements from the proper perspective since it is tacitly 
assumed that the matter content of the model does not differ from the minimal one. This 
is not so in the true string model which possesses several fields beyond the minimal ones 
at poorly determined intermediate mass scales. Nevertheless, no well motivated choice for 
these mass scales has improved the situation, and there even exists a no-go theorem to 
this effect (see below). 

Two ways out this impasse have been proposed ]3"B |: string threshold effects on the 



gauge couplings could reduce the effective unification scale down to C(10 16 GeV) []34|,[35 
or new matter fields could delay unification until O(10 18 GeV) [|32| , |36f| . The first alternative 
is not viable in models built in the free fermionic formulation since threshold corrections 



always increase the effective unification scale [pl| , p4| .n The minimal field-theoretical ex- 
tra matter content (beyond the fields in the supersymmetric standard model) needed in 
the second alternative has been determined to be an extra pair of Q and D c vector-like 



representations |3~T j, which fit inside a 10 representation of SU(5). However, in a more 



complicated situation like for a string- derived flipped SU(5) model, one can only assert 
that one extra Q, Q pair is absolutely necessary. The main purpose of this paper is to 
search for such models within the chosen formulation. 

The organization of this paper is as follows. In Sec. 2 we review the fundamentals 
of the free fermionic formulation and rephrase its model-building rules to suit our compu- 
tational purposes. In Sec. 3 we describe our computer algorithms. In Sec. 4 we explore 
the parameter space for the desired models, and in Sec. 5 we present the results of the 
extensive computer search. In Sees. 6 and 7 we describe in some detail the most promising 
model found, and in Sec. 8 we summarize our conclusions. 



2 The threshold corrections needed in models where this mechanism may work require rather 



unnatural values of the moduli fields [34]. 



2. The free fermionic formulation 

It has become well known that consistent classical string vacua (i.e., "string mod- 
els") in low spacetime dimensions (d < d c = 10), especially in four dimensions, can be 
constructed in a variety of formulations |l|. Different formulations correspond to dif- 
ferent choices for the two-dimensional world-sheet conformal field theory, so that these 
extra two-dimensional degrees of freedom, which may or may not admit a spacetime inter- 
pretation, together with the two-dimensional conformal fields that describe the observed 
four-dimensional spacetime, ensure the (super) conformal invariance and modular invari- 
ance of the full string theory. In the free fermionic formulation |||3|,^|,[5],|6|,[F] , all the inter- 
nal two-dimensional degrees of freedom are fermionized by utilizing only free world-sheet 
fermionic fields. The stringy consistency conditions are cast into a set of constraints on 
the spin-structures |3^] of all the world-sheet fermions (i.e., on the boundary conditions of 



these fermions as they are parallell transported around noncontractible loops on the world- 
sheet) and a set of constraints on the relative contributions of different spin-structures to 
the string partition function. 

The two-dimensional world-sheet fermions should provide the exact amount of central 
charge for the Virasoro algebra in order to render the quantum theory free of conformal 
anomalies, both in the left- and right-moving sectors. For the four-dimensional heterotic 
stringi, cl = 9 and c# = 22 are needed; these conditions fix the total number of world- 
sheet fermions. In the light-cone gauge, in addition to the two left-moving fermionic fields 
ijj^ (fi = 1, 2) which are the superpartners of the two transverse string coordinates X^, the 
internal fermionic content consists of 18 left-moving and 44 right-moving real fermions, 
respectively. The world-sheet supersymmetry in the left-moving sector is nonlinearly real- 
ized among the left-moving fermions which must transform as the adjoint representation 



of a semi-simple Lie algebra of dimension 3(10 — d) = 18 [j39 j. To achieve N = 1 space- 
time supersymmetry, this algebra must be SU(2) 6 [|7|], and therefore we can group the 
left-moving fermions into 6 triplets (% , y l , u ) (£ = 1, 2, . . . , 6), each transforming as the 
adjoint representation of SU(2), with the following supercurrent 

6 

T F =^d 2 X M + ]TxW, (2.1) 

i=i 



3 We follow the convention that left-movers are supersymmetric and right-movers are non- 
supersymmetric. 
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where Tp must be periodic or antiperiodic on the world-sheet, corresponding to spacetime 
fermions or spacetime bosons. 

In addition to the requirement of (super) conformal invariance, the most important 
constraints on the spin-structures of world-sheet fermions are obtained by considering 
modular invariance. The basic idea is to investigate the properties of the string partition 
and correlation functions under modular transformations, and then search for the solutions 
that accommodate modular invariance. In the free fermionic formulation, besides the con- 
straints that modular invariance imposes on the one-loop string partition function, further 
constraints can be obtained in two different approaches: (1) requiring modular invariance 
and factorization of multiloop partition functions 0,0]; or (2) performing a "physical sen- 
sible projection" from the space of world-sheet states onto the subspace of physical states 
which ensures the correct spin-statistics connection These two approaches have been 
shown to be entirely equivalent |)]||. It turns out that the results of such analysis can be 
described as a set of model-building rules. We now give a self-contained account of these 
rules. 

The first set of rules restricts the spin-structure assignments of the world-sheet 
fermions. Before stating these rules, it is helpful to introduce the spin-structure vectors 
which specify unambiguously the spin-structures of all world-sheet fermions. Consider a 
(/-loop world-sheet S g ; there are 2g noncontractible loops. For each loop, a spin-structure 
of a set of world-sheet fermions is given by a orthogonal matrix representation of the 
first homotopy group 7r(£ 3 ) which is non-abelian for g > 2. In the approach of Ref. H, 
since only the one-loop partition function is relevant, one can just specify the boundary 
conditions of the world-sheet fermions around the two noncontractible loops of the torus, 
and they are mutually commuting. However, in the approach of Refs. [0,|3| , on a two- or 
higher- loop world-sheet, the analysis of multi-loop modular invariance and factorization 
becomes a rather delicate affair, since non-commuting boundary conditions are allowed || . 
Nevertheless, simplification can be achieved by restricting oneself to the subset of mutually 
commuting boundary conditions, which amounts to replacing 7r(E g ) by its abelianized ver- 
sion, i.e., the first homology group i?i(S 5 ). Under this restriction, even in the approach of 
Refs. 0,|3| ( which we follow closely in this paper) the spin-structure matrices can be simul- 
taneously diagonalized for all noncontractible loops. As a result, for any noncontractible 
loop v, one can simply specify the corresponding spin-structure by a spin-structure vector 
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where n = n r + n c , n r and n c are the numbers of the real and complex fermions. The 
component v(f) represents the boundary condition of fermion / under parallel-transport 
around this loop 

/ ^ _ e *™(/)/. (2.3) 

It is convenient to introduce a sign for each spin-structure vector v as follows, 5 V = e inv (^' J '), 
This sign distinguishes between spacetime bosons (5 V = +1) and spacetime fermions (5 V = 
— 1) [0|§. Clearly, one can choose to define the spin-structure vector differently. In fact, 
the W-vector used in Ref. ||] and the F-vector used in Ref. are related to our f-vector 
byv(f) = l-2W(f) = -2V(f). 

It was shown in Ref. || that if one only considers the special spin-structure vectors 
whose components are rational numbers, then all possible rational spin-structure vectors 
that give rise to consistent string models must form a finite additive group H, which can be 
generated by a basis B = {&i, 62, • • • , b p } of a finite number of spin-structure basis vectors. 
Each basis vector bi has order N bi , which is defined as the smallest positive integer such 
that 

N bi bi = (mod 2). (2.4) 

Here is a special spin-structure vector which assigns antiperiodic boundary conditions 
to all world-sheet fermions; we will refer to it as the Neveu-Schwarz vector in the rest of 
this paper. The basis B must obey the following constraints: 
(Al) It can be chosen to be a canonical basis, such that 

p 

^771^ = (mod 2)^^ mi = (mod N b .) VI (2.5) 
i=i 

In terms of this canonical basis, an arbitrary spin-structure vector v G S is an integer 
linear combination of basis vectors reduced in interval (—1, 1], which we express as 



i=l 



0,l,2,...,(JV 6i -l). (2.6) 



(A2) One of the basis vectors, which we choose to be the first one 61, satisfies 

^N bl b t = l (mod 2). (2.7) 

Here 1 is another special spin-structure vector which assigns periodic boundary conditions 
to all world-sheet fermions; we will refer to it as the Ramond vector for convenience. 
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Although the general forms of such b\ vector have been tabulated we will choose 
bi = 1 for simplicity, hence N± = 2. 

(A3) If N bib . is the least common multiple of N b . and N b . : then 

N b . b .(bi-bj) =0 (mod 4). (2.8) 

(A4) For any bi e B, 

( (mod 4) if N b . is odd; 

N bi (b i -b i ) = \ 1 (2.9) 

I (mod 8) if N bi is even. 

However, if N bi {bi + 1) = (mod 4), then bi must satisfy a stronger condition 

N b .(bi-bi) = N b .(l-1) (mod 16). (2.10) 
In (|2.8|)-( pTT0| ), the dot product between two spin-structure vectors u and v is defined 

byi 

w - v = Ue+ E ~lJ2- E }«(/)«(/)■ (2- 11 ) 

^ real complex real complex ' 

left left right right 

(A5) The number of real fermions that are simultaneously periodic in any given three 
basis vectors bi,bj,bk is even ||, i.e., 

EW)WM/)=0 (mod 2) (l<i,j,k<p). (2.12) 

real 

(A6) The spin-structure corresponding to each bi should be an automorphism of the Lie 
algebra defining the world-sheet supercurrent, and all such automorphisms must commute 
with one another. In the case of our interest, the world-sheet Lie algebra is SU(2) 6 , and 
the possible commuting sets of automorphisms have been tabulated in Ref. [0. This rule 
restricts in a rather intricate way the possible sets of basis vectors which are allowed to 
coexist. We defer to the next section a further discussion of this constraint. 

Having reviewed the constraints of the spin-structure vectors, we now move on to the 
second set of rules which restricts the relative contributions of different spin-structures to 
the one-loop string partition function. As it was shown in Refs. [|3]||, this set of rules can 
be given in terms of some proper normalized coefficients which enter the one-loop partition 

4 This differs from the definition in Refs. by a overall sign due to the different conventions 
used. 
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function and depend upon the spin-structures. Since such one-loop spin-structure coeffi- 
cients are pure phases ||, any direct computer implementation would involve complex 
numbers. In order to avoid this unnecessary complication, we use a /c-matrix notation sim- 
ilar to that introduced in Refs. ||||. More precisely, for a given basis B = {&i, 62, • • • , b p }, 
we define a p x p positive integer matrix (which we call the /c-matrix) in terms of the one- 
loop spin-structure coefficients C[ 6 *] associated with basis vectors 6j and bj (see Ref. @), 
such that its element kij is given byi 

= 6 ii e 2 * ik "/ N "i, (l<fc y <N b .). (2.13) 

With the /c-matrix, the second set of model-building rules can be written as follows: 
(Bl) There is freedom in choosing k\\. In fact, we have 1 < k\\ < N bl and 

8/cii = N bl (h -h) (mod4iV bl ), (2.14) 

or 

8fcn = N bl [h ■ h - 4&i(^)] (mod 8JV 6l ). (2.15) 

Since we take b\ = 1, Eq. ( |2.14| ) gives kn = 1 and Eq. ( |2 . 1 5| ) gives k\x = 2. 

(B2) Besides kn, all other diagonal elements ka (i > 2) of the /c-matrix are completely 
fixed by the corresponding kn, 

8fc« = JV 6i [4fe il + 6 i -6 i + l-l-2(2-JV 6l )6 i (^)] (mod 8iV 6i ) (i > 2). (2.16) 

Note that the last term on the right side vanishes when N bl = Ni = 2. 

(B3) The elements below the diagonal kij (i > j) and the corresponding elements 
above the diagonal are related as follows (for i 7^ j), 

A(N bi hj + N hj k 3l ) + 2N h N bj [biW) + bj(^)] ~ N bi N bj (b t ■ b 3 ) = (mod 4N bi N b .). 

(2.17) 

It is easy to see that a /c-matrix is completely determined by kn and the elements 
below the diagonal (or equivalently the elements above the diagonal). Therefore, for a 
given basis B, there are 2]^^ - g bibj distinct /c-matrices associated with it [[|, where g bibj 
(i 7^ j) is the greatest common divisor of N bi and N b . (which gives the number of choices 
for kij). 

5 Note: Despite its same name, our /c-matrix definition differs from that in Refs. [f|||. 
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Given a pair (B, k), subject to rules (Al)-(A6) and rules (B1)-(B3) respectively, one 
can construct a consistent four-dimensional heterotic string model. The total Hilbert space 
of the model is a direct sum of the sub-Hilbert spaces of all the sectors in the model, each 
of which corresponds to a spin-structure vector v of form (|2.6|) and is refered to as the 
v-sector and sometimes denoted as H v . Furthermore, the physical states in the v -sector 
are determined by the generalized GSO-projections [3D], which take various forms |3|,|J. 



Strictly speaking, in order to perform the GSO-projections correctly, special care must 
be taken when dealing with real fermions, because the zero modes of real fermions with 
periodic boundary conditions contribute to the string partition function ||. Of course, 
if all the left- and right-moving real world-sheet fermions in the model can be separately 
complexified (i.e., described by left- and right-moving complex fermions obtained by pair- 
ing up two real fermions that have the same spin-structure) then the model is just a special 
case of a complex fermion model, and there is no real subtlety involved. 

In this paper, however, we would like to consider a more general class of models, 
in which only 2ul (2til < 20) of the left- and 2n# (2n# < 44) of the right-moving real 
fermions are complexifiable in the usual sense. Clearly, the two left-moving fermionic fields 
should always be treated as one complex fermion. (Our choice of ^-vector (described 
below) allows complexification of the x into three complex fermions H.) In addition, 
any remaining left- and right- moving real fermions appear in equal numbers, and can be 
grouped in left-right pairs which share the same boundary conditions (periodic or antiperi- 
odic) for all basis vectors, effectively representing a non-chiral conformal field theory, e.g., 
an Ising model || . This class of pseudo-complex fermion models is a subset of the general 
class of real fermion models considered in Ref. || , and is precisely the class of "real fermion 
models" discussed in Ref. [|J . Restricting ourselves to the pseudo-complex fermion models, 
we can divide the world-sheet fermions into three categories: (1) left-moving complex 
fermions; (2) ur right-moving complex fermions; (3) nj Ising models each formed from a 
"residual" left- and a "residual" right-moving real fermion. The numbers of these three 
types of world-sheet fermions are related by 

m = 44 - 2n R = 20 - 2n L . (2.18) 

Accordingly, we can rearrange any spin-structure vector v into three parts: a n^- 
dimensional vector v L , a n#-dimensional vector v R and a n/-dimensional vector v 1 . For the 
particular class of basis vectors we will consider in this paper, below we will see that rule 
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(A6) implies that the elements in v L and v 1 can only be either 1 or 0, whereas the entries 
in v R can be rational numbers, which describe genuine right-moving complex world-sheet 
fermions. For any complex world-sheet fermion /, one can define the fermionic charge in 
the f-sector as [IH^Tl 

Q v (f)=F v (f) + lv(f), (2.19) 

where F v (f) is the fermion number and v(f) the boundary condition. Therefore, for 
left- and right-moving complex fermions, there are the fermionic charge vectors Q% (Tri- 
dimensional) and Q R (n^-dimensional) respectively. For the real fermions in left-right 
pairs, the concept of the fermionic charge is not applicable. However, in the models under 
consideration, their effect can be taken into account simply by introducing two operators 
for each Ising model / : (1) a "fermion number difference" operator Ai ? „(/ / ), which equals 
the fermion number of the left-moving real fermion minus that of the corresponding right- 
moving fermion when v{f I ) — 0, but equals zero when v(f ) = 1; (2) a "chirality" operator 
r^(/ J ) such that T„(/ 7 ) = 1 if vif 1 ) = 0, and r„(/ r ) = ±1 if v{f J ) = 1. Collectively, we 
have two well-defined nj-dimensional vectors AFj and T^. 

In terms of Q L , Q R , AF 1 and T 1 , we can rewrite the GSO-projections in a very 
convenient way. In fact, for all basis vectors hi (1 < i < p), the physical states in the 
f-sector must satisfy 

b^Q^-b R -Q R +b(-AF^ + ^bi-(Ti-l I )=v(r)+h(r)+2j2^ ( mod2 )> (2- 2 °) 

where l 1 is the n/-dimensional vector with all the entries equal 1. When a particular string 
state of f-sector that satisfies (|2.20|) is obtained, its quantum numbers (with which the 
corresponding vertex operator can be constructed) are completely determined by the three 
vectors Q%, Q R and §. More explicitly, the first entry of describes the spacetime 
spin of the state, the next three entries give the charges under the three £7(1) 's which 
make up the conserved Uj(l) current of the N = 2 world-sheet supersymmetry algebra, 
and the remaining entries (til —4) are the charges under residual left-moving global U(l) 
symmetries. The right-moving fermion charge vector Q R encodes partially i the charges 

6 It is quite often the case that a full matter representation of the gauge group is formed by 
states that come from more than one sector of the model. 



10 



under the Kac-Moody currents and therefore the gauge quantum numbers. In particular, 
the rank of the total gauge group G equals the dimension of Q^, that is 

rank(G) = n R . (2.21) 

Finally, the vector is related to the eigenvalues of the order (disorder) operator for all 
the Ising models, which we denoted collectively by a n/-dimensional vector al with 

a v (f I ) = v(f I )T v (f 1 )- (2-22) 
The mass of a string state in the f-sector is given by the following formula || 

M 2 = -i + + W/M/) + E( fre ^ uencies )' 

^ real complex * left 

left left 

= - 1 + l{lJ2+ E }^(/)^(/) + E( frec i uencies )- ( 2 - 23 ) 

^> real rnmnlnx * vin-bt 



real complex ' right 

right right 



Where the oscillator frequencies for a fermion / are 

l + v(f) 



+ integer, (2.24) 



2 

and if / is complex, its complex conjugate /* has frequencies 

1 - + integer. (2.25) 
Also, the bosonic oscillators contribute integer frequencies. 

3. Description of the computer program 

One of the advantages of the free fermionic formulation is that the rules for model- 
building are quite simple and can be computerized, which makes a fairly systematical study 
of the string models built in this formulation possible. We are aware of one such attempt 
in the literature [[IT], which was based on the formulation of Ref. [|J. However, in Ref. |TT 



only complex fermion models were considered, while the more interesting models with real 
fermions were not taken into account. In order to investigate the pseudo-complex fermion 
models, we have independently developed a computer program, based on the formalism 
described in Sec. 2 and with an approach different from that of Ref. JITJ . Our computer 



software is written in FORTRAN. The bulk of the program is devoted to the generation of 
the massless spectrum of a pseudo-complex fermion model once the basis and the consistent 
fc-matrix of the model are given. In addition, our program is also able to calculate the 
complete trilinear superpotential and identify all the fermion Yukawa couplings of the 
model. In this section we explain our method and discuss some algorithms involved. 
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3.1. The basis and the k-matrix 

To obtain a basis B that obeys the rules (Al)-(A6), we should first decide on the 
first basis vector b\. An obvious choice that we make throughout this paper is b\ = 1. 
As long as b± is fixed, rule (A2) is not relevant, and the remaining task is to generate 
additional basis vectors such that the other rules are satisfied. We choose the second basis 
vector 62 to be the "supersymmetry generator" S, which allows the existence of spacetime 
supersymmetry by giving rise to massless gravitinos. Of the several possible forms of S 
that are compatible with the choice b\ = 1 0, we choose the simplest one, which is the 
only one that has been used for mo del- building in the literature, namely 

S = (1 100 100 100 100 100 100 : R ). (3.1) 

Here the first entry is reserved for the two transverse ip^ treated as a single complex fermion; 
the following 18 entries correspond to the left-moving internal fermions treated here as 6 
triplets (x^, c</); and the vector Or after the colon, which separates the left-movers 
from the right-movers, simply indicates that all right-moving fermions have antiperiodic 
boundary conditions. 

An immediate bonus of choice (|3.1| ) is that the possible forms for the left-moving part 
of all other basis vectors can be rather easily determined. In fact, if S is given as in ( |3.1|) , 
then in order for the basis vectors to satisfy rule (A6), the left-moving internal fermions 
must have either periodic or antiperiodic boundary conditions |7[], organized as 6 triplets 
(x 1 V ) ui ) °f rea l fermions. Furthermore, in any basis vector bi, the boundary conditions 
for these fermions must satisfy [[7] 

W) = Ux") + h{y l ) + h{J) (mod 2), (3.2) 
biiWMtfMv^Mu') ^{0,1} (1<*<6). 

Because of our choice of S, the otherwise rather intricate check of rule (A6) can be done 
very easily in the class of models we consider. In addition, the specific form of the vector 
S also gives us some useful information on the right- moving part of all other basis vectors. 
From (|3.1| ) it is easy to see that only y £ and u (£ = 1, ... ,6) are available to serve as 
the "residual" left-moving real fermions (since the x £ can always be complexified). In 
order to obtain the pseudo-complex fermion models, we need to allow equal numbers of 
left- and right-moving real fermions. Therefore, we deduce that the first 12 entries of the 
right-moving part (corresponding to 12 real fermions denoted by y £ and u e ) of all basis 
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vectors will be either or 1. Furthermore, we can always treat the remaining right- moving 
fermions as complex, i.e., the last 16 entries of each basis vector can be rational numbers 
in the interval (—1,1]. 

Any test vector with the above general features has a chance of being a viable basis 
vector. Once such a vector is generated (either by hand or by computer) our program 
first checks whether or not it satisfies rules (A3) and (A4), by calculating the dot product 
of this vector with all the other basis vectors that have been previously entered into the 
basis. If the test vector fails this step, the program discards it and repeats the same 
procedure for another test vector, until a good candidate is found. At this point, the 
program records the relevant information about this new basis vector, such as its order 
and the associated dot products. This whole process can be repeated as many times as 
one wishes and is controled by an interactive command. After having obtained a set of 
several basis vectors, this set would constitute a basis B if rule (Al) and (A5) are also 
satisfied. The algorithm for checking rule (Al) is very simple. Our program first generates 
all possible linear combinations of form ( |2.6| ) with the set of basis vectors, and then looks for 
the Neveu-Schwarz vector 0, which should appear only once, when all the "coordinates" 
Vi = 0, if the basis is canonical. As far as rule (A5) is concerned, it can simply be 
checked using (|2.12| ). However, since we are particularly interested in the pseudo-complex 
fermion models, we have designated a special subprogram to check whether the model is 
complexifiable as described in Sec. 2. Since all complexifiable models satisfy rule (A5), 
once this is established we would have obtained an allowed basis B which would give rise 
to a pseudo-complex fermion model. It is worthwhile to point out that in our approach we 
start with basis vectors of the general form described above, and do not specify ul, ur or 
nj. These numbers are determined by the program after the complexification procedure 
is successful. The advantage of this approach is that our program can treat equally any 
pseudo-complex fermion model with gauge group of rank between 16 and 22 (see Eqs. 
( [Hp and (gUID). 

Having generated a basis B, all consistent /c-matrices can be systematically obtained 
by using ( |2.13| )-( |2.17| ). In practice, the number of possibilities is prohibitively large and 
one has to resort to random generation of these matrices (see Sec. 5). However, even 
in the random approach there are still many distinct /c-matrices, although they do not 
necessarily lead to distinct or phenomenologically acceptable string models. To avoid a 
large cosmological constant at the string scale, we are interested only in models with 
spacetime supersymmetry. Furthemore, to obtain chiral fermions this better be N = 
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1 supersymmetry (as opposed to N = 2 or 4 which are also possible). Such models 
can be obtained if the basis B contains enough basis vectors to break the initial N = 4 
supersymmetry down to N = 1 However, since the GSO-projections relevant to the 
massless gravitinos in the model (we need one and only one gravitino for N = 1) depend 
on the choice of /c-matrix, only a subset of the /c-matrices will give N = 1 supersymmetric 
models in the end. We have written a special subprogram which calculates the massless 
gravitinos provided by the S'-sector immediately after a k- matrix is generated. In this 
way, the program acts on only those selected /c-matrices and discards all others. In Sec. 
4.2 we will show that given a specific basis, the constraints on the /c-matrices imposed by 
requiring N = 1 supersymmetry can be worked out explicitly, in which case we just use 
the explicit constraints to test for N = 1 supersymmetry. 

3.2. The massless sectors 

To find the massless spectrum of a string model, we only need to consider those sectors 
which could possibly contain massless states. All we need are the basis B and the mass 
formula ( |2.23| ). For an arbitrary v -sector, let us consider the following two quantities 

v L ■ v L + -v 1 ■ v 1 , (3.3a) 
V R. V R + v i. v i_ (3.35) 

right right 

In the models under consideration, since the left-moving part of all basis vectors is made 
up of O's and l's, the left-moving part of any spin-structure vector v of form (|2.6| ) will only 
have or 1 entries as well. Therefore, there are only two possible forms of massless states in 
the left-moving part: the non-degenerate bosonic vacuum with a single fermionic oscillator 
of frequency | acting on it, and the degenerate fermionic vacuum with no oscillator at all. 
Accordingly, a i>-sector will provide massless states only if Ml = or Ml = 4. 

For the right-moving part, since the last 16 entries can be arbitrary rational numbers 
in the interval (—1, 1], a similar analysis is to certain extent basis-dependent. Nevertheless, 
if we confine ourselves to models in which the only rational numbers that appear in the 
basis vectors are ±|, then these would also be the only rational numbers in any vector v. 
As a result, the possible right- moving fermionic frequencies are 4, \ and |, and in order 
for a f-sector of such models to provide massless states, it has to have Mr = 0, 2, 4, 6, 8. 



^ real complex * 

left left 

^ real complex * 
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3.3. The gauge group 

The massless vector states of a heterotic string model, which are created by exciting the 
left-moving antiperiodic transverse fermions ^ and the right-moving internal fermions, are 
the gauge bosons. @ These come from massless sectors with Ml = and transform as the 
adjoint representation of the gauge group G of the string model. Since the right-moving 
complex fermions are the building blocks of the two-dimensional Kac-Moody currents 
that underlie the four-dimensional gauge symmetry, the gauge group G of each pseudo- 
complex fermion model can be systematically and unambiguously determined by studying 
the fermionic charge vectors Q R for all the sectors with Ml = 0. First of all (as mentioned 
in Sec. 2), the number of right-moving complex fermions ur, which only depends on the 
basis £>, gives the rank of G. Secondly but most importantly, the Q R vectors which label 
the gauge bosons always constitute the roots of the Lie algebra associated with the gauge 
group G §. 

We start with the gauge bosons in the Neveu-Schwarz sector (0-sector) which is always 
present in all models (see (|2.6|)). Due to the existence of the "residual" antiperiodic right- 
moving real fermions (not present in complex fermion models), to get gauge bosons we 
can either excite two right-moving complex fermions or excite one right-moving complex 
fermion and one such "residual" right-moving real fermion (in both cases we have Mr = 0, 
i.e., '00' gauge bosons). In the first case we get right-moving fermionic charge vectors of 
the form 

Q R (f i )=±S ij ±S ik (j^k), (3.4) 
as well as some vanishing ones. In the second case, these vectors take the form 

Q R (f i ) = ±8 ij . (3.5) 

The Neveu-Schwarz sector is special in that the Q R vectors of the gauge bosons in this 
sector themselves constitute the roots of a Lie algebra associated with a subgroup of G 
with of same rank, which we refer to as the "minimal" gauge group G m of the string 
model. Clearly, for the simple models in which the Neveu-Schwarz sector is the only sector 
giving gauge bosons, G m and G coincide. It is remarkable to see that, from ([T3]) and 

7 Gauge bosons which arise by exciting left-moving internal degrees of freedom preclude the 
existence of chiral (massless) spacetime fermions if the resulting gauge group is abelian (non- 
abelian) ||. 
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dT5|), not only the Q R vectors correspond to the roots of G m , but their elements also 
specify precisely the coordinates of these roots in the standard orthonormal basis fl4~3 |. In 
addition, the forms of ( |3.4| ) and ( |3.5| ) imply that the minimal gauge group G m in general 
is a direct product of simple factors of SO (2m), SU(m + 1) and SO(2m + 1), as well as 
some U(l) factors. Here, the non-simply-laced groups SO(2m + l) arise simply because the 
"residual" right-moving real fermions are allowed in pseudo-complex fermion models (see 
( |3.5|) ). Since roots always appear in pairs of opposite sign, we only need to consider the 
positive roots, i.e., the Q R whose first nonvanishing entry is positive. In the Neveu-Schwarz 
sector, our program first constructs the 462 positive roots of the Lie algebra D22, which 
corresponds to the primordial gauge symmetry SO(44) of four-dimensional heterotic string 
models built in the free fermionic formulation. The positive roots of the minimal gauge 
group Gm are then obtained by performing the GSO-projections ( |2.20|) which in this case 
reduce to 

-bf ■ Qq +bf ■ AFq = (mod 2). (3.6) 

Because the Qq vectors are njj-dimensional, and in the pseudo-complex fermion models 
we have 16 < < 22, only the first ur entries of the 462 positive roots of D22 enter in 
( p.6| ). (When tir < 22, the first ur entries of some of these positive roots are all zeros; 
such n^-dimensional zero vectors are obviously not positive, and our program discards 
them automatically.) After the positive roots in the Neveu-Schwarz sector are found this 
way, the program is instructed by an interactive command, to either find the simple roots 
from these positive roots (according to the algorithm which we will describe below) and 
therefore pin down the minimal gauge group G m for some testing purposes, or to go on 
searching for all the additional gauge bosons in other sectors and finally find the total 
gauge group G of the model. 

Besides the Neveu-Schwarz sector, any f-sectors of the model with Ml = could 
also provide gauge bosons. The most common gauge bosons of this type have Mr = 8 
(i.e., '08' gauge bosons); in practice we have also found models with '06' gauge bosons. 
The procedure to find these additional gauge bosons is the same: we start from a set 
of appropriate fermionic charge vectors Q R in each of such sectors, keeping only those 
vectors whose first nonvanishing entry is positive, and then find the subset which survives 
the GSO-projections. In contrast with the Neveu-Schwarz sector, two differences are worth 
mentioning. First, for the Neveu-Schwarz sector the GSO-projections ( |3.6|) are independent 
of the /c-matrix. So once a basis B is given, the minimal gauge group G m is completely 
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determined (prior to the generation of the /c-matrices). For other sectors that give gauge 
bosons, the corresponding GSO-projections become 

-bf .Qf + 6j.AF/ = 2^^ (mod 2), (3.7) 

3=1 bj 

which clearly depend on the /c-matrix. As a result, the total gauge group G of the string 
model varies with different choices for the /c-matrix. Second, although the gauge bosons 
of the Neveu-Schwarz sector give rise to the minimal gauge group G m , the gauge bosons 
from other sectors by themselves do not have this property. Instead, the additional gauge 
bosons transform as certain representations of the minimal group G m , such that all these 
representations of G m combine with the adjoint representation coming from the Neveu- 
Schwarz sector to make up precisely the adjoint representation of an enlarged group G. 

Having found all the positive roots, to uniquely identify the gauge group G, the next 
step is to find the simple roots and calculate the Cartan matrix. The same idea was followed 



in Ref. ||41|1 , where a recursive procedure was used to find the simple roots according to 
the property of their inner products with each other. Here we take a simpler approach, 
which is based on the following definition of the simple roots: A simple root is a positive 



root that cannot be written as the sum of two other positive roots [43]. Our program first 
adds the positive roots together pairwise, then it compares all the positive roots with these 
two-root pairs and finds the simple roots. This way we end up with n srt < ur simple roots 
in some arbitrary order. 

To determine explicitly the factors in the gauge group we proceed in two steps. We 
first isolate a set of the simple roots which forms a subspace by itself, i.e., that would give 
a block diagonal submatrix in the Cartan matrix in a suitable basis. That is, ai -ctj ^ for 
i,j G b, but «i ■ ctfc = for i e b, k G" b. The number of roots in this block is the rank 
of the sought after gauge group factor Gi. We then calculate the Cartan matrix Cij of this 
block of simple roots. If r % = 1 and C X1 = (2) then G l = U(l) (SU(2)). For r % > 1 the 
determination of Gi can be made unambigously from and det C, as seen from Table I, 
and the following additional considerations: for simplicity we identify Ci with Bi, although 
proper identification is possible; we distinguish B-j from E-? by the fact that the Cartan 
matrix of the former is not symmetric; and we take into account explictly the isomorphism 
A3 ps D 3 . The above procedure is repeated until all simple roots are exhausted and then a 
symbolic expression for the total gauge group is printed out. In order to be able to identify 
the matter representations automatically, it is necessary to reorder the simple roots in each 
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block in a standard way (such as the one given in Ref. |[45| ). This is possible once the Gi 
have been identified. However, the algorithm is non-trivial and cumbersome and will not 
be described here. 

When n srt < ur, the gauge group G also contains (ur — n STt ) U(l) factors. In the 
models that we are considering, the U(l) gauge group can arise in two ways: (1) If a 
single right-moving complex fermion has unique boundary conditions so that no other 
right-moving complex fermions have the same ones in all the basis vectors, then this single 
complex fermion provides a U(l) factor, and the corresponding U(l) charge of a given 
representation is specified by the value of its fermionic charge in that sector; (2) If a unitary 
group arises through the embedding SU(m) x £7(1) C SO (2m), then a U(l) appears as 
a "trace", and the U(l) charge of a given representation in this case is the sum of the 
fermionic charges over all fermions in the trace. However, there are more complicated 
embeddings for U(l) factors, identification of which needs to be done on a non-automated 
case-by-case basis. 

3.4- The matter states 

The massless matter states can come from massless sectors with Ml = or 4. Similar 
to the case of massless gauge bosons, the gauge group properties of the massless matter 
states are carried by their right-moving fermionic charge vectors Q^, which are the weights 
of the corresponding representations of the Lie algebra associated with the gauge group 
G Q. (As explained in Sec. 2, the left-moving fermionic charge vectors also contain 
useful quantum numbers of the states ||.) 

Since we only consider models with N = 1 spacetime supersymmetry, the matter 
states form N = 1 supermultiplets. A full N = 1 supermultiplet is provided by states from 
a f-sector and the corresponding (v + S')-sector, so we always treat them together. The 
fermionic and scalar components of an A = 1 supermultiplet transform under the same 
representation of the gauge group, but their left-moving quantum numbers are different. 
For the fermionic components we choose to keep states with positive spacetime helicity, 
i.e., Qv(^) = \- The correct left -moving fermionic charges for the corresponding scalar 
components then satisfyi 

£ Qv(fi)Qgrav(fi) = i (3-8) 
i=2,3,4 

8 This result is a generalization of Eq. (37) in Ref. || [with the following notational rela- 
tions: Qv(fi) «>/3,7; Qgrav(fi) *-* ai2,a34,a56] for the case of an arbitrary component of the 
gravitino vertex operator. 
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where Qg rav is the fermionic charge vector of the (positive-helicity) gravitino in the model, 
which our program has determined, and the sum runs over the three complexified left- 
moving fermions constructed from the x l '■ 

The way in which we find the right-moving fermionic charge vectors for all mass- 
less matter states is essentially the same. We first construct a set of appropriate fermionic 
charges vectors for each sector, and only keep those that survive the GSO-projections 
( p.20|) . Since the weights of Lie algebra representations do not necessarily come in pairs, 
we have to consider both positive and negative weights. To specify the irreducible rep- 
resentations of the matter states, we only need to keep track of the highest weights that 
appear in the end. Again, differently from the approach of Ref. ||41|1 , we find the high- 
est weights by first converting all weights into the corresponding Dynkin labels using the 
properly reordered simple roots, and then keeping only those weights whose Dynkin labels 
are all non- negative integers | ]45|| . 

The identification of the transformation properties of the selected set of weights under 
the various Gi groups is done via a look-up table. That is, the program knows the Dynkin 
labels of all allowed representations under all gauge groups which can occur in a level-one 
Kac-Moody construction (a limited set |10[)- Looping over all Gi we obtain a symbolic 
label for each state, e.g., 10-1-1 for a 10 under SU{h) and a singlet under other two GVs. 
For later purposes, the multiplicity and conjugacy classes [f45l , |TC)|1 (see Sec. 3.5) of each 



representation are also calculated and saved. We also determine the contribution to the 
anomaly of each SU(N), N > 3 subgroup Gi by each of the states in the spectrum, as 
follows 

A l = J2 ai(Rj) dimiRj/Gi), (3.9) 
R j 

where a(R) = (N—3)l(N—2n)/(N—n—l)\(n—l)l for totally antisymmetric representations 
(the only ones that can occur at level one |l(|]) of dimension dim(R) = ( ) |46]. In Eq. 



dim(Rj/Gi) is the dimension of Rj under the subgroups other than Gi. Since all 
these anomalies must vanish, this constitutes a rather valuable check of the correctness of 
the program. 

3.5. The superpotential 

Once the full massless spectrum is available, the program computes the cubic super- 
potential couplings using the techniques of Ref. [§]. This is done by forming all pos- 
sible </>i</>2</>3 sets and determining whether such coupling is invariant under the various 
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(3.10) 



two-dimensional world-sheet symmetries. One of these is equivalent to four-dimensional 
gauge invariance. This is checked using the conjugacy classes determined earlier: un- 
der Gi we must have Ci(<pi) + c^cfo) + q(03) = Omod(dj). For example, for SU(n), 
c(l) = 0,c(n) = l,c(n(n- l)/2) = 2, c(R) = -c(R), and d = n; for SO(2n), 
c(l) = 0, c(n) = 2,c(2 n_1 ) = —1, and d = 4. We also compute explicitly the confor- 
mal field theory correlators for the various Ising models (corresponding to each left-right 
real fermion pair). For cubic couplings these can be tabulated explicitly using the following 
multiplication table 

I a + cr_ / 

I f I a + a- f \ 

a+ a+ I f cr_ 

a- a- f I a+ 

f \ f a- °+ I J 

where a± are the eigenvalues of the chirality operator defined in Eq. ( |2.22| ) (also known 
as the order/disorder Ising model operators), and / represents the corresponding left- or 
right-moving real fermion. For example, the correlator (a+a-f) is nonvanishing because 
a + o~-f according to Eq. (|3.1Q|) . 

For flipped £77(5) models, among the various cubic couplings there are some of the 
form: 10 a -5-5, 10*, ■ 10 c ■ 5^, 5 ■ 1 ■ 5 e . These are candidates for the up- and down-quark 
type and lepton-type fermion Yukawa couplings. The code identifies these couplings and 
determines whether it is possible to have non-vanishing top, bottom, and tau Yukawa 
couplings at the cubic level, i.e., whether there are couplings such that 10 a = 10& or 
10 a = 10 c , and 5 d = 5 e . 

A separate portion of the code determines the set of possibly non-vanishing quartic 
and higher-order superpotential couplings according to the (significantly more complicated) 
rules of Ref. . All we want to point out here is that the multiplication table given above 
can be used to determine whether a high-point Ising model correlator (corresponding to 
real fermion degrees of freedom not involved in the picture changing operations) vanishes. 
This test reduces considerably the number of terms that need to scanned by hand. In fact, 
at the quartic level this program gives only non- vanishing couplings. 
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4. The space of flipped SU(5) models 

The computer program described in the previous section can be used to study generic 
models built within the free fermionic formulation. As discussed in the Introduction, we 
restrict ourselves here to models with the flipped SU (5) group as the observable part of 
the total gauge group. This constraint is rather weak once one realizes the size of the 
space to be explored. Instead of doing an exhaustive search of all possible bases and k- 
matrices, we choose a well motivated subset of the space of basis vectors (seven vectors) and 
consider a general class of vectors 'a' which effect the breaking SO(10) — > SU(5) x U(l), 
as discussed below. We also consider all the possible /c-matrices compatible with the bases 
under consideration. 



4-1. The choice of basis 

We study a basis of eight vectors B = {1, S, bi, 62, &3, &4, 65, a}, the first seven of which 
are given in Table II. The choice of 1 and S was discussed in Sees. 2 and 3.1. The presence 
of vectors b\ and 62 is necessary to reduce the spacetime supersymmetry from N = 4 
down to N = 1. The precise form of these vectors is basically fixed by self- consistency 
constraints []T|,|4"T|j. The presence of 63 is necessary to enlarge the massless spectrum and 



possibly accommodate three generations of chiral fermions. It has even been argued |)47 



that the bi, bi-, 63 set is a unique one which must be present in any three-generation model. 
Basis vectors 64 and 65 are needed to allow for symmetry-breaking Higgs representations 
in the massless spectrum. Their precise forms have not been fully classified, although their 
left- moving entries are constrained by self-consistency conditions (see Sec. 2). Also, these 
vectors should not destroy or enhance the spacetime supersymmetry of the model. In sum, 
given the form of S, the first five vectors in B are to a large extent "written in stone", i 



whereas the forms of 64 and 65 have been chosen as in the revamped model |2C 



The vector a serves the purpose of breaking the resulting S'O(IO) gauge symmetry 
down to SU(5) x U(l), and must therefore contain rational entries which we take for 
simplicity to be \. We study the following generic form of the a- vector 

a=(00yV ■•• Oj/V : ffffffw^wW^w* \\\\\ \\\ ||||H00), (4.1) 

where Eq. (|3.2| ) has been used to set a(w e ) = at(y ). The choice for the fixed entries 
was made to a large extent to follow the revamped example. With the choice of 61,2,3,4,5, 



The other possible forms of S @ have remained virtually unexplored in the literature. 
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the complexification procedure indicates that (depending on a) the possible gauge groups 
will have rank 16-19. The actual rank in excess of 16 depends on how many of the three 
£7(1) symmetries left unbroken by 61,2,3,4,5 are also left unbroken by a. These symmetries 
are generated by the following complexified pairs of real fermions, U a : w 2 w 3 , U(, : y 1 w 6 , 
and U c : y 4 y 5 . It is convenient to divide up the possible forms of a into eight classes, 
according to the number and type of U aj b, c symmetries that they break or do not break. 
This classification is given in Table III. (The revamped model belongs to class 3a.) 

4-2. The constraints on the k-matrix 

Given the basis B with eight vectors, the corresponding /c-matrix is 8 x 8, but only 
the 28 lower-diagonal entries are independent (see Sec. 2). As shown below, this leads 
to a very large number of possible k- matrices, most of which may not necessarily lead to 
phenomeno logically interesting models. We know present three classes of constraints on 
the k- matrix which guarantee N = 1 spacetime super symmetric models, and allow the 
possibility of three generations of matter fields and symmetry-breaking Higgs fields. 

The GSO projections in Eqs. (|2.20|) can be simplified considerably for potential grav- 



itino states in the S'-sector, as follows 

bf-Q^ = l + bi(^) + k i2 (mod 2). (4.2) 

For bi = a we have a L ■ Qg = (since a(x e ) = 0) and a(^^) = 0, therefore /c§2 = 1- This 
constraint follows solely from the assumed form of a in Eq. Ql.l|) . Analogously one can 
show that b\ ■ Qg = b\ ■ Qg and b% ■ Q$ = 65 ■ Qg, and therefore /c62 = &32 and kj2 = ^42 if 
we want to preserve at least one gravitino. It can also be shown that /C52 is determined as 
follows /C52 = 1 + &21 + &32 + &42 (mod2). If we pick the positive-helicity gravitino state, 
i.e., Qsi'tp' 1 ) = o) then the eight possibilities for Q^x 1 ' 2 , X 3 ' 4 , X 5 ' 6 ) follow uniquely from 
the values of £32,42,52 as follows: Q£(x 12 ' 34 ' 56 ) = - |(-l) fc32 - 42 ' 52 . 

We must also have representations left from the 61,2,3-sectors. In this case a special 
property of these vectors, namely bi ■ bj = ^(tp^) — Ylk=i M^) = — 4 (where tp^ 1 are the 
first five right-moving complex fermions) allows one to deduce that 

5 

bk ■ Q£ - b* ■ Qg = 6f • Q£ - b* ■ Qg = Q£ (r) ~ £ Og (^)- (4-3) 

k=i 

It can then be shown that we must have /C43 = £53 for the states in the 61-sector to remain. 
A similar exercise for the 62,3-sectors yields the constraint /C43 = /C53 = £54. Further study 
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of the effect of the remaining GSO projections indicates that if /C43 = hs = /C54 = 1 (2), 
then we will get 16,16' (16,16') representations from each of 61,2,3- 

An analogous argument applied to the 64^-sectors yields further constraints on the 
fc-matrix, as follows 

hi = k 6 3 hi = h 5 , hi 7 1 h3 => hi ^ &65, (4.4a) 
k 71 = k 74: k 73 = k 75 , k 7 i ^ k 74i => k 73 ^ k 75 . (4.4b) 

4-3. The possible models 

Given our choice of basis B = {1, S, 61, b%, 63, 64, 65, a}, it is convenient to determine 
the possible sets of SU (5) representations coming from the sectors Ti-bx^MM,^ smce these 
will contribute the Ramond states of the observable spectrum. It is not hard to show that 



the possible representations are as follows: 

61,62,63 : 5 + 10 or 2 -10 or 2 -5, (4.5a) 

64, 65 : 5 + 10 or 5 + TO or 10 + TO or 5 + 5, (4.56) 

64 + 65 : (5 + 5)„ or 2 • 5 or 2 ■ 5, (4.5c) 



where the 5,5,10,10 representations are pieces of 16,16 representations of SO (10), whereas 
(5 + 5)„ = 10 of SO(10). We can now form all possible combinations of the above represen- 
tations such that the total SU (5) anomaly is zero,0 and nio > 4, rij^ > 1. Five possibilities 
arise: 7110^^:4/1,4/2,5/1,5/2,6/2. (These go with n 5 /n§ : 2/5, 2/4, 1/5, 1/4, 0/4.) Of 
these possibilities, only 4/1 and 5/2 models are likely to yield three light 10's (or equiva- 
lently n g = n\o — rij^ = 3). The important point here is that the 4/1 and 5/2 models are 
obtained only when the (5 + 5)„ representations are the ones contributed by 7^64+65- This 
requirement imposes a constraint on the possible a- vectors as follows. 

Since the 5O(10) symmetry is not broken down to £77(5) x U(l) until the a-vector 
GSO projection is performed, whether we get (5 + 5)„ or 2 • 5 or 2 • 5 from the 64 + 65- 
sector depends only on this last GSO projection. Reverting to the old notation for the GSO 
projection, the coefficient which matters is C[ b4 ~^ b5 } = ±1 (±i) for (64+65) -a = even (odd). 
By explicit calculation one can show that if C[ 64+65 ] = ±i, then (5 + 5)„ will survive if 
a (y 2 ) = a (y 3 ) (y 2 y 2 an d y 3 y 3 form two Ising models), otherwise we get 2 • 5 or 2 • 5. Given 

10 Otherwise a non-trivial mixture with the hidden sector must occur. (States from the NS- 
sector contribute zero to the anomaly.) 
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the form of a in Eq. ( f4.1| ), it follows that a(w 2 ) = a(w 3 ) and therefore w 2 w 3 can be 
complexified and the U a : w 2 w 3 symmetry remains unbroken. Analogously, one can show 
that if C[ &4 ^ &5 ] = ±1, then (5 + 5),; survives if a(y 2 ) ^ oi(y 3 ) and therefore a(w 2 ) 7^ a(w 3 ) 
and U a is broken. In sum, if (64 + 65) ■ a = even (odd) then U a will be broken (unbroken). 
Now, 

(b 4 + b 5 )-a = -l + i{a(y 2 )-a(y 2 )+a(y 3 )-a(y 3 )+a(w 1 )-a(w 1 ) + a(ty 4 )-a(tZ; 4 )} = -1, 

(4.6) 

since y 2 y 2 , y 3 y 3 , w 1 w 1 , w^w 4 always form Ising models. Therefore U a must remain unbro- 
ken and from Table III we see that only a- vectors from classes l,2a,2b, and 3a will respect 
this. The universe of possibilities then consists of a-vectors from these classes and the 
set of distinct /c-matrices compatible with each basis. Furthermore, only 5/2 models will 
contain the extra pair of Q, Q representations required for gauge coupling unification. 

5. The computerized search 

Given the large number of possibilities, it is clear that one will benefit from a statis- 
tical approach to model generation. However, even random generation of a statistically 
significant number of /c-matrices (for a fixed a) is prohibitive: a basis of p vectors goes with 
anpxp /c-matrix, whose p(p — l)/2 lower diagonal entries are independent. In our case 
p = 8 and there are 2 28 « 268 million possible /c-matrices. Fortunately the constraints 
on the /c-matrix derived in Sec. 4.2 are very powerful. Extensive random generation of 
/c-matrices indicates that there are further redundancy factors to be obtained, as follows: 

(i) N = 1 supersymmetry determines ^52,62,72,82 once /c2i,32,42 are given. However, the 
eight choices for ^21,32,42 (which correspond to the eight possible choices for the su- 
persymmetry generator) give equivalent models. We then get a 2 7 reduction factor. 

(ii) To obtain states from 7~£b 1 ,b 2 ,b 3 we need to impose /C43 = /C53 = /C54 = 1, 2. If the model 
ends up being a 5/2 (or 2/5) model, the specific choice just picks 5/2 over 2/5 (i.e., 
interchanges 10 <-> 10). This gives a 2 3 reduction factor. 

(iii) There is no obvious analytic constraint on the entries ^31,41, 51, and heuristically we 
find that their values do not affect the phenomenological properties of the models.0 
We therefore fix them; a 2 3 reduction factor. 

11 They may even give rise to equivalent models, as expected from the internal discrete sym- 
metries associated with the various Ising models. 
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(iv) Finally, we have the analytical correlations among k^i^^^s and £71,73,74,75 to obtain 
states from TLb^bs- A reduction factor of 2 2 . 

The combined analytical and heuristical reduction factor is 2 15 which reduces the 
universe of possible /c-matrices down to 2 13 = 8192. This number of fc-matrices is in 
practice still quite large and makes random generation of /c-matrices belonging to this 
class practically necessary to be able to explore several ct-vectors. A word of caution: since 
the number of possibilities is finite, one must make sure that each generated fc-matrix 
is distinct from all previous ones.0 This is not hard to do; we simply identify each k- 
matrix by the decimal equivalent of the 13-digit binary number formed by its independent 
entries in a fixed order. In order to gain confidence on statistical results (e.g., percentages 
of models with a given feature) obtained by sampling ^ 10% of the possible models, we 
generated all 8192 /c-matrices for one a-vector and found the exact results to be in complete 
agreement with those obtained statistically from smaller samples.lll 

The a-vectors which we consider introduce a large dimension in the space of possible 
models. Each of these a- vectors has 24 free real fermion entries (2 24 possibilities), although 
these are highly constrained. Rule (A6) (see Eq. ( |3.2|) ) requires ct(y £ ) = a(w ) for this 
type of a- vectors (they all have ct(x e ) = 0). There are also six more constraints on the 
24 entries, one from rule (A4) and five from rule (A3) (the constraint from bi = 1 follows 
from rule (A4), and that from bi = S is automatically satisfied for our choice of a- vector). 
Finally, the splitting of a- vectors into classes (see Sec. 4.1) introduces correlations among 
three pairs of entries. Therefore, there are 24-6-1-5-3=9 independent real entries in any 
given a-vector, a total of 2 9 = 512 possibilities. 

The program described in Sec. 3 is capable of producing a large amount of output. 
It was therefore necessary to screen possible models, i.e., calculate the whole spectrum 
and cubic Yukawa couplings, and keep only summary information about each model (e.g., 
fc-matrix number, gauge group, mo /tijq, fermion Yukawa couplings, number and type of 
hidden sector fields, etc.). With the k- matrix information it was then easy to retrieve any 
particular model afterwards. 

12 This is not very important for small samples of fc-matrices since these obey the following 
probability distribution n dlstinct = N(l - e ~ nt " als/N ) sa n trials + 0(n 2 trials /N) (N = 8192 in this 
case) . 

13 This run took ~ 40 hours of cputime in a Silicon Graphics 340S computer. In comparison, 
the 2 28 A;-matrices would take « 153 years of cputime! 
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The results of our search with the general class of a-vectors we consider can be sum- 
marized as follows: 

(a) Class 3a. There are two rank-17 gauge groups possible for any a-vector in this class 
(the fc-matrix determines which one occurs): SU(5) x SU(4) x SU(2) 4: x U(l) 6 and 
SU(5) x SU(4) x SO (10) x U(l) 5 . We did not find any 5/2 model with the a-vectors 
explored and for all possible choices of /c-matrices. The revamped model (a 4/1 model) 
is a prototype of models belonging to this class. 

(b) Classes 2a and 2b. There are two rank-18 gauge groups possible for each a-vector: 



SU(b) x SU(4) x SU(2) 4 x U(l) 7 and SU(b) x SU(4) x 50(10) x U(lf. All 5/2 



models found (these occur 1/8 of the time; about 1000 per a- vector choice) come with 
the second gauge group. They also have the same number of massless fields (67) , and 
the same number and type of hidden matter fields (six 4,4 and seven 6 of SU(4), and 
three 10 of SO (10)). These models fall into two subgroups: half of them have a 2/3/2 
Yukawa set (2 up-quark-type, 3 down-quark-type, and 2 lepton-type cubic Yukawa 
couplings) and the other half have a 1/3/2 Yukawa set. 
(c) Class 1. There is a large number of rank-19 gauge groups possible in this case, 
although they do not occur for all a- vectors (some are more 'prolific' than others). 
The following is a list of the 13 gauge groups which we have been able to identify: 



It is interesting to note that (contrary to naive expectations) the addition of the a- 
vector to effect the breaking down to flipped SU(5), can quite likely enlarge the final 
gauge group and leave no such group after all. There are two kinds of a-vectors in 
this class which give 5/2 models: 




SU(9) 



x 



SO(10) x U(lf 
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(cl) The 'price'-like a-vectors, one of which was used in Ref. []32| to construct the first 
5/2 flipped SU(5) model in the literature. In this case 5/2 models (which occur 
1/4 of the time) always come with the gauge group SU(5) x SU(5) x SO(10) x 
C7(l) 6 and do not possess right-handed leptons in their massless spectra (they 
also have problems with the anomalous £^4(1) cancellation), and are therefore 
irrelevant for phenomenological purposes 

(c2) The other type of a- vectors give 5/2 models (1/8 of the time) with gauge group 
SU(5) x SU(4) x SO(10) x U(l) 7 and have no problems with right-handed leptons 
or cancellation of £/a(1)- These models have the same hidden matter spectrum 
as those 5/2 models belonging to classes 2a and 2b, and they also split into two 
subgroups: half of them have 1/3/2 and the other half have 1/4/2 Yukawa sets. Ill 



6. The phenomenology of the various models 

In the previous section we exhibited the list of possible 5/2 models found under our 
assumptions about the basis vectors. Given the Yukawa couplings of the various models 
one can easily convince oneself that only the models that belong to classes 2a or 2b which 
have the 2/3/2 Yukawa set are of any interest. Indeed, in models with only 1 up-quark-type 
Yukawa coupling (i.e., 1/3/2 or 1/4/2) one must identify that coupling (F t fh) with the one 
for the top quark, since no other such coupling appears at cubic, quartic, or quintic order 
in the superpotential. The problem is that Ft does not appear in terms of the form FfFh 
at cubic, quartic, or quintic order, and therefore we do not expect a (sizeable) bottom 
Yukawa coupling. We thus concentrate on the 2/3/2 models (in comparison, the fairly 
successful revamped model is a 2/3/3 (although 4/1) model). 

14 In Ref. [ p2| | the gauge group for the model presented was incorrectly determined to be 
SU(5) x SU(A) x 50(10) x C7(l) 7 , and so were some of its massless representations. This problem 
may have arisen because contributions of additional gauge bosons belonging to new '08' sectors 
may have been overlooked. 

15 There is a further subdivision inside each of these sets into two different kinds of 1/3/2 and 
1/4/2 Yukawa sets. 
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6.1. The model 



The model we consider in some detail here belongs to class 2b (although all 2/3/2 
models in classes 2a and 2b are equivalent to this one) and has the following a-vector 

a = (00 000 000 000 000 000 011 : 000001 011001 \\\\\ \\\ Ji±±1100), (6.1) 

and /c-matrix (this is the same one as the one for the revamped model, although there are 
512 such fc-matrices that will give the same results) 

/2 1 2 2 2 2 



1111 

2 2 2 2 
2 2 2 
2 2 2 
2 2 2 
2 2 2 



1 1 

2 2 



2 
2 
2 
2 



2 
1 
2 
2 
2 
2 
2 
2 



4 
1 
1 
4 
1 
3 

3/ 



(6.2) 



\2 1 1 1 1 1 

We start by listing in Table IV the massless fields and their transformation properties 
under the rank 18 gauge group SU(h) x U(l)~ x SO {10) h x SU(4) h x U(lf. The cubic 
and quartic superpotential terms are easily calculated ||. We obtain 

W 3 = gV2{F F 1 h 1 + F 2 F 2 h 2 + F 4 F 4 h x + F 4 f 5 h 45 + F 3 f 3 h 3 + f 2 l c 2 h 2 + f 5 l c 5 h 2 

+ ^4^503 + \F A F&o + F A Fji x + F 5 F 5 h 2 

+ (hi~h 2 ®12 + h 2 h 3 $ 23 + h 3 hi$ 31 + h 3 h 45 (j) 45 + h.c.) 
+ ^(045045 + + + + 0~0~ + <Pi<Pi + h 45 h 45 )§ 3 + (VlV2 + VlV2)$0 
+ (^304 + 03</>4)$5 + ($12*23*31 + $120 + 0~ + $120*0* + h.C.) 
+ ^$31+^3 $31 

+ AsA>$23 + L»iL> 2 $23 + D 5 D 5 $ 23 + D 7 D 7 $ 31 + D 3 D 3 ^ 31 

+ ±D 5 A3$0 + ^D 5 D 7 fa 

+ F 4 F 6 $ 12 + \F 3 f 4 <S> + \F 2 F^ 3 + F 6 F 4 (f)+ + 75^404 
+ F x F 2 D b + F 2 F 4 r 2 }, 



(6.3a) 



and 

W 4 =F 2 f 2 h 45 4 + F 3 F 4 D 4 D Q + F 3 F 5 D 4 D 7 

+ l c 3 F 3 F 6 D 7 + ltF 2 F 3 4> 3 + F 1 l , 3(0 + 03 + 0"0s) 
+ P 3 F 5 D 7 0- + F 2 F 5 D 3( p- + F 2 F G D 3( p 4 + F h F x D 2 D 7 
+ F 5 F 2 DiD 7 + F 3 F 3 D 3 D 6 + F 4 F 3 D 4 D 7 + F 5 F 4 D 5 D 7 . 
Calculable coefficients M have been omitted from W 4 . 



(6.36) 
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6.2. F- and D-flatness 

To preserve unbroken supergravity at the string scale, the D-terms of all U(l) sym- 
metries and all the F-terms must vanish, i.e., 

dW 

(W) = < — > = 0, (6.4a) 

0(pi 

<Ai>=X>iK&>| 2 + e = 0, (6-46) 

i 

(D a ) = J>*| (^)l 2 = 0, a = 1, 2, 3, 4 , (6.4c) 

i 
i 

where e = g 2 Tr £/a/1927t 2 . In these formulas, q A is the charge of the i-th field under the 
anomalous U A (l) g8|: U A = (l/TrU A ) ELiP^] U U 9a are the charges under the 5-1=4 
orthogonal traceless linear combinations, and Yi are the charges under (the anomaly-free) 
U(l)~. We obtain 

XJ a = (-36C7i - 1217a + 24U 3 - 12[/ 5 )/Tr ^, Tr U A = 46.4758, (6.5) 



and thereforeE^I e = (3.7 x g x 10 17 GeV) 2 . Note that Tr t/4 = 0; in fact, almost all singlets 
(except 771,2,771,2) are neutral under this ?7(1) symmetry. The rotated traceless Z7 (l)'s 
(excluding U4) can be written as follows (proper normalization is understood, although 
not relevant in what follows) 

Ua = -3U 1 -U 2 + 2U 3 -U 5 , (6.6a) 
U{ = U 3 + 2U 5 , (6.66) 
\j' % = U x - 3U 2 , (6.6c) 
U£ = 3Ut + U 2 + 4U 3 - 2U 5 . (6.6d) 

The constraints on the singlet vevs obtained by solving Eqs. Q6.46, c| ) are as follows 

x 45 = e/15-^ 3 2 , (6-7a) 

x+ -x- = e/15+ \v£, (6.76) 
x 3 \ - x 2 3 = e/5, (6.7c) 
^ (a*) + 2(x 23 + x+ - X12) = V\~ V 2 2 , (6.7a 7 ) 

i=3,4 
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To recover the proper mass units we recall that we have set k = \/&it/Mpi = 1 in Eq. ( |6.46 ). 
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where x p = |(0 P )| 2 — |(0 P )| 2 , and Vi = (z^ c ), Vi = (z^ c ). Similarly, .D-flatness of the 
anomaly- free U±(l) and U(l)~ requires 



|Vbr = |Vir + 2(|r?ir - l^il + \m\ - \m\), (6-8a) 

i=0 



= ^|Vi| 2 = |1/ 4 | 2 + \V 5 \ 2 = V 2 . (6.76) 



The F-flatness constraints (Eqs. ( |6.4a| )) can also be worked out for the cubic super- 
potential and the resulting 24 equations (one for each singlet field) can be summarized as 
follows 

$12*31 = $12*23 = 045 $3 = and h.c, (6.9a) 
$23$3i + + 0~ + 4>i4>i = and h.c, (6.96) 

045045 + 0i0i + + + + 0~0~ = 0, (6.9c) 

0304 + 0304 = 0, (6.9d) 
+ $3 + 4>~<ti2 = 0~$3 + + $i2 = and h.c, (6.9e) 
5¥ 4 + TO+M2 = 0, (6.9/) 

-^V 4 V 5 + 3 $3 + 04$5 + 20 3 $12 = 03*3 + 04*5 + 20 3 $12 = 0, (6.9o) 

4 $ 3 + 3 $ 5 + 20 4 $i2 = and h.c, (6.9/i) 
07i,*7i,»72,*/2}$o = 0. (6.9i) 

These equations have five classes of solutions when solved in conjunction with the D- 
flatness constraints in Eqs. (|6.7| ). There are eight possible solutions to Eqs. ( |6.9a|) , i.e., 



$12 = $12 =$3=0, $12 = $31 = $23 = $3 = 0, $12 = $12 = 045 = 045 = 0, • • •. Four of 

these solutions violate either of the D-flatness conditions in Eqs. (|6.76, cj) . The remaining 
five solutions can be worked out in detail and are as follows: 

1. $12 = $12 = $3 = 0, 

a. $ 5 = 0, F 4 = 0, (|3), (|6.96, c,d, /I) ; 

b. $ 5 = 0, V 4 ± 0,F 5 = 0, (|3), (16.96, c, d, /I) ; 

c. $ 5 + 0, V 4 = 0, 03,4 = 03,4 = 0, Q, flS.96, c,/l ); 



d. $ 5 ^ 0, V A ^ 0, 03,4 = 03 = 0, O, ( P-96, C,/,fl| ) 



2. $12 = $31 = $23 = $3 = 0, 0+ = 0" = , 03,4 = , 045 045 = 0, V 4 = 0, v\ > V% , 

a. $ 5 = 0, (U, (S3); 

b. $ 5 ± 0, 03,4 = 0, (|3), (PQ) . 
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3. $12 = $31 = $23 = $3 = 0, 0+ = <f>~ = 0, 3i 4 = 0, 045<^45 = 0, V\ < 

a. $5 = 0, (§3), ( |B.96,/,flD ; 

b. $ 5 + 0, 03 = 0, Q, (16.96, Q - 

4. $ 12 = $31 = $23 = 0, 045 = ^45 = 0, 0+ = 0" = 0, V? = 2e/15, 

a. $ 5 = 0, 03,4 = 03,4 = 0, V\ > V 2 2 , O, Q6.96,/,g| ); 

b. $ 5 ^ 0, (13), ( |5.9&,c,d, Z"j"J ). 

5. $12 = $31 = $23 = , 045 = <^45 = 0, 0+ = 0" = 0, F 3 2 = 2e/15, 

a. $ 5 = 0, 03,4 = 04 = 0, V\ < Vl ([TTj), ( |B.96,/,flO ; 

b. $ 5 ^ 0, flB.96, c,d, Z"j"J ). 

In cases 2 and 3, the solution to 045 045 = depends on the sign of £ = e/15 — ^V 3 2 . For 
£ > (< 0), 045 = (045 = 0), whereas for £ = , 045 = 045 = 0. Note also that cases 
2 and 3 are mutually exclusive. In solving Eqs. (|6.9z| ), we have assumed that 771,2,^1,2 



have generically nonvanishing vevs and therefore $0 = 0. If 771,2 = 771,2 = 0, further 
simplification of the above solutions is possible. 

6.3. Symmetry breaking and EVA mechanism 

There are five high-energy scales in this problem, which are correlated and in principle 
self-consistently determined, as follows: 

(i) the string unification scale M v « 7.3 x g x 10 17 GeV JEJggfDI 5 

(ii) the scale of XJ a anomaly cancellation vevs (0) ~ = C(10 17 GeV); 

(iii) the SU(5) x U{1)~ -> SU(3) C x SU(2) L x U(l) Y symmetry breaking scale M x , 
undetermined theoretically but phenomenologically of O(10 15 ~ 18 GeV); 

(iv) the scale where the surplus Q,Q get heavy M eva (the EVA scale H50[]), to be deter- 
mined; 

(v) the scales of hidden sector condensation A4 and A10, to be determined from the hidden 
sector matter spectrum. 

The gauge symmetry breaking occurs as follows 

SU(5) x U(l)~ x *7 4 (1) d SU(3) C x SU(2) L x U{1) Y > x U{1)~ x *7 4 (1) 

-> SU(3) C x Stf(2) L x U(1) Y , (6.10) 

where we have assumed that U^l) remains unbroken, i.e., 7712 = 771,2 = 0. This ap- 
pears to be a reasonable choice since these fields do not participate in the anomalous 
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Ua(1) cancellation. Also, if 771,2,771,2 7^ 0, the subsequent discussion will remain qualita- 
tively unaffected. We will generally assume that Vq,Vi,V5 7^ 0, V4 = 0, and optionally 
^2,3, 7^ 0. Given the charges of Vq 1? under the relevant £7(1) symmetries, it is clear 
that only Y = \V - \Y remains unbroken below the scale Mx- In fact, there are two 
flat directions of the scalar potential: cpi = \{i>q + vf) + ^75^5 an d ^2 = ^75(^0 + v x)- The 
imaginary parts of the scalar components of these fields are eaten by the broken U(l) gauge 
bosons and the remaining supermultiplet components become heavy higgs/higgsino states. 
There is another linearly independent field which remains light, the so-called flaton/flatino 
supermultiplet [^TJ. The <ig x ,d% components of Fq,i,Fs appear in the higgs triplet mass 
matrix (see Subsec. 6.4). We are then left with the Qo,i,Qa,5 components. The scalars 
get either eaten by the X, Y gauge bosons or become heavy higgs bosons, whereas the 
fermions interact with the X, Y gauginos through the following mass matrix 







Q 5 


Y 


Qo 




(5) 




M 1/2 = Qi 


(4) 
W[ 


(5) 
W[ ' 


95V1 


X 




95V 5 


/ 



(6.11) 



where come from higher-order superpotential couplings, and effect the EVA mech- 

anism. We expect all these fields to become heavy, although at different mass scales (see 
Sec. 7). 

6.4- The Higgs mass matrices 



The Higgs doublet masses originate from the hihj 
couplings and to quartic order are 



HiHj and Fifjh k 



ViLjH k 
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v 4 J 



(6.12) 
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The Higgs triplet mass matrix which effects the doublet / triplet splitting receives contribu- 



tions from: h{hj —> 
resulting matrix is 



DiDj, FiFjhk 



V % D k d% FiFjhu 



VidfDk, and FiFj 



d)d\. 



The 







D 2 


D x 
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(5) 
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<f>3 ) 



(6.13) 



6.5. A possible numerical scenario 



Until the full set of quintic (and possibly higher-order) superpotential couplings is 
calculated, one cannot work out some features of the model such as the precise eigenvalues 
of the doublet and triplet mass matrices, or the diagrams contributing to proton decay, 
or the details of gauge coupling unification in this model. Nevertheless, we can postulate 
a possible scenario (to be confirmed by further calculations) given the limited amount of 
information available at this time. 

There are three basic phenomenological considerations which help decide on the vevs 
left undetermined by the flatness conditions: (i) dimension-five operators in proton decay, 
(ii) the fermion Yukawa couplings and the higgs doublet mass matrix, and (iii) the higgs 
triplet mass matrix. Cubic couplings of the form F a Fbh D Q a QbD, F c f h D Q C LD, 
hh<p D DDcj) would be disastrous for proton decay [^6[]. We have (see Eq. ( |6.3a| )): 
F Q Fxhx,F 2 F 2 h2,F 4 F 4 hx; F 4 f 5 h 45 , F 3 f 3 h 3 ; hxh 3 ^3i, h 2 h 3 ^ 23 . With our assumptions in 
Sec. 6.3, Qo,x do not contain quark fields, whereas Q2,3,4 do. It is then important to have 
$23 = 0, and perhaps also $31 = 0, assuming Q 2 (Q 4 ) contains second (third) generation 
quark fields. It turns out that these two constraints are either automatic or can be imposed 
(and still maintain flatness) in the five F-flatness cases in Sec. 6.2. In case 1 both Hx and 
H 2 become pure massless states (see Eq. ( |6.12| )), whereas in cases 2 and 4 (3 and 5) Hx 
(H 2 ) is a pure massless state. Based on the sensible assumption that (so far neglected) 
higher-order contributions to Ai 2 will likely make massive (light) mixed states rather than 
(light) pure states (since the former have more possible couplings), we conclude that cases 
2 and 4 are the most promising ones, since the others would give too many cubic fermion 
Yukawa couplings. 
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Let us analyze the two preferred cases in turn: 

(a) Case 2. We need to impose $23 = 0, while we get $31 = automatically. To avoid 
H45 becoming heavy, we need the 045 = solution to 045045 = (which assumes 
that £ = e/15 — ^V 3 2 > 0, and is then consistent with the choice V 3 = 0) and V 2 = 

2 — 

(also consistent with V 5 > V 2 ). We then get Hi,H 2U oc -($> 31 )H 2 + ($> 12 )H 3 , 
£2,3,5; H 2 , #45 hght and H 23h oc ($i 2 )H 2 + ($ 3 i)#3, H 45 ; H u H 3 heavy. We expect 
higher-order contributions to M. 2 to give an intermediate scale mass to one pair of 
the remaining light doublets. The cubic fermion Yukawa couplings which remain 
are then -F4/5/145 D Qau^H^ and F^F^hi D Q^d^Hi, leading to the identification 
At = Xb = gV%- At the quartic level we have (see Eq. ( |6.36|) ) cF 2 f 2 h^(f)4 which 
will give a vanishing coupling since 04 = in this case. Preliminary probing into the 
quintic couplings reveals a structure of the form 

F 2 F 2 /i 1 {0+0-, 0303, (6.14a) 

MM0+0", 0303, 0404}, (6.146) 
/ 5 /5^l{^303,0404}. (6.14c) 

All these couplings will also vanish since + = 0~ = 03 5 4 = in this case. This 
scenario appears disfavored. 

(b) Case 4- Here we also need to impose $23 = 0, whereas $31 = is automatic. To 
obtain a light H45 we need to set $3 = V 2 = V4 = 0. However, V 3 = 2e/15 > in this 
case and the higgs doublet mass matrix has a novel structure. We get Hi, H^, H 23 £, 
L 2 ,5', H 2 ,H^ 5 light, and H 23 h,L 3 ; Hi,H 3 heavy. Note that there will be a mixing 
between "higgs" and "lepton" doublets in this case, since we expect one pair of the 
remaining light doublets to become heavy when higher-order contributions to M. 2 are 
included. The advantage of this case over Case 2 above is that the vevs contributing 
to the higher-order couplings are generally non-zero. We obtain: At = Xb = gV%, 



\ c = C (0 4 )/M, with c = O(l) and (<f>)/M ^ 1/10 PI- Of the three quintic couplings 
in Eq. (|6.14|) , the first of these could give the strange-quark Yukawa coupling, the 
second the tau Yukawa, and the third one the \i Yukawa. For this assignment to be 
realistic we would need the second coupling (A T ) to be sizeable despite its potential 
suppression by (0) /M, whereas the first coupling should remain small (and so should 
the third one). These details may be possible to arrange since there are three sources 
of uncertainty at this point: the actual coefficient of the terms, the size of the singlet 
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vevs involved, and the possibility of mixing between Q2 (d?,) and some of the other 
Q's (d c 's) in the symmetry breaking process (higgs triplet mass matrix). 

A study of the triplet mass matrix M.3 is not very illuminating at this time, although 
it is clear that we must demand V4 = $0 = </>3 = to obtain a light d%. It is also clear that 
both tig x become heavy. In case 4, D\2] Di, d§, and a linear combination of D2, -D3, d\ are 
also heavy. Higher-order contributions to .M3 are needed to continue this analysis further. 



7. Gauge coupling unification 

Let us now address the question which motivated our extensive search for a 5/2 model. 
If one starts from the low-energy values of the gauge couplings (sin 2 6 W = 0.2331 ± 0.0013 
|H , as = 0.113 ± 0.004 p3| , and a" 1 = 127.9 ± 0.2) and expects to get gauge coupling 



unification at a scale Mjj, then various intermediate-scale particles have to contribute in 
suitable ways to the running of the gauge couplings. The SU(5) x U(l)~ breaking scale 
Mx is an unknown in the problem, but not the only one since (as discussed in Sec. 6.3) 
there are other scales inter- related to it. As an example we could imagine that Mx is close 
to Mjj , in which case the expressions of Refs. [[32],|37|] apply (for Mx < Mjj see Ref. [p6|), 
as follows 

L Q = (23.31 - 25.86) + \L V . + \L E , + 2L~, (7.1a) 
L D c = (58.29 - 67.90) + L L + \L E , + 2L~ - \L Uc , (7.16) 

where Lr = Y2i ln(Mj/ /to_rJ and the sum runs over all the supermultiplets in representa- 
tion Ri. The constant ranges allow for values of the low-energy parameters inside the 1-cr 
error ellipsoid p7 |. In the present model there are no U c or E c representations besides 



the standard ones, thus L\j c = Le c = in ( |7.1|) . The important point encoded in these 



expressions is the fact that Lq > and therefore intermediate-scale extra vector-like Q 
representations are needed. These are not present in 4/1 models, and one pair is present 
in 5/2 models. Note also that L^c > (since Lu c = 0) and we must also have extra 
vector-like D c representations at intermediate mass scales, but this is always the case in 
flipped SU (5) models. 

From Table IV and the vev choices made in Sec. 6, the contribution to the running 
of the gauge couplings can be split up into five groups: 
1- -^2,3,4, ^2,3,5, ^2,3,5 contribute the usual three generations of quarks and leptons. 
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/2 \ ('2') 

2- Hi, 2,3,45) -Hi.2,3,45 contribute to = ln(M;y /m^ ) 2 , where are the eigenval- 
ues of ■ (Note that includes the contribution from the expected two light 
higgs doublets.) 

3. D ll2l3) 45,< 5 ; ^1,2,3,45,4,1 contribute to L D < = EX^/m{ 8 >) 2 , where ^ are 
the eigenvalues of .M3-M^ with the last three columns removed from 

4. Qo, Qi^Qh = (Qo + Qi)/V2, Qeva = (Qo - Qi)/V2 and Q 4 , Q 5 ^Q H = (F 4 Q 4 + 
Vr Q^)/V,Q eva = (V^Q A - V±Q 5 )/V, where Qh,Qh become heavy higgs/higgsino 
states in the 577(5) x U(l)~ symmetry breaking process (see Eq. ( |6.11| ) with w <C 
V,V) and Qeva, Qeva get masses M eva = 0(w). (Note that in Eq. ( |6.11| ) we need 

7^ or 7^ to avoid a massless Qeva, Qeva pair.) The latter contribute 
to LQ = 2ln(Mu/M eva ). 

5. Fi, Fj which have ±| electric charges contribute to L~ = Y2i hi(M;y/m^ ±1/ ' 2 ' ) ) 2 where 

(±!/ 2 ) are e ig envames f ^ ne mass matrix with FiM.~F 

From Eq. ( pMa|) it is clear that Lq > (23.31 — 25.86) since L~ > (recall that Ljjc = 
Le c =0). From the expression for Lq just derived, it follows that M eva & 10 12 GeV. This 
upper bound is reduced to « 10 10 ' 8 ' 6 GeV for M x = 10 17 ' 16 ' 15 GeV |fj. This is the only 
direct test of gauge coupling unification in this model which we can perform at this time 
since it only involves the determination of the scale M eva . We have explored the quintic 
terms in the superpotential and found terms of the form ^0,1-^4,5(0) (DD , FF) /M which 
would give w ~ M eva ~ O(10 5 GeV) if the SU(4) h condensation scale is C(10 12 GeV) and 
{(f)) ~ 10 17 GeV. These are plausible numbers which do not violate the upper bound on 
Meva even for Mx = 10 15 GeV, and still leave room for a non-vanishing L~. 



8. Conclusions 

Gauge coupling unification is one of the few universal predictions of string models 
(together with the presence of gravity and gauge interactions). An even more pervasive 
fact is that the string unification scale can be precisely calculated in any given string model. 
The robustness of this prediction is in sharp constrast with the basically unlimited number 
of possible models with or without supersymmetry, and with all possible gauge groups and 
matter representations. (The space of models possesses certain internal structure which 
for example correlates gauge groups and matter representations.) Nevertheless, string 
unification in any given model can actually be tested. Indeed, any such model gives in 
principle definite predictions for some of the best measured parameters at the electroweak 
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scale, namely sin 2 9 W and a^. From this point of view, traditional unified models are also 



testable this way ||33|| , although there the unification scale has to be chosen to fit the 
experimental data. 

The motivation for this paper was to search for a flipped SU(5) model which could 
possibly accommodate the measured values of the low-energy gauge couplings or equiv- 
alently which could unify at the string scale. There is a necessary condition for such a 
model (in level-one Kac-Moody constructions and with small string threshold effects), that 
it possesses extra Q, Q representations, i.e., a 5/2 model. To pursue this objective within 
the free fermionic formulation, we have developed a sophisticated and comprehensive com- 
puter program which can sweep systematically over large numbers of models and determine 
their gauge group, massless spectrum, and cubic superpotential, and decide to keep only 
"interesting" models. Our search space, even though vast, is by no means exhaustive and 
much of it still remains uncharted. Nevertheless, the structure of what we did explore has 
proved to be very simple. Of the few supersymmetric 5/2 models which we found, only 
one of them allows for both top- and bottom-quark Yukawa couplings. Further probe into 
the superpotential and the D- and F-flatness constraints indicates that the model is phe- 
nomenologically sound. However, much work remains to be done to bring this model up to 
the level of development of its predecessor (the revamped model). An encouraging result 
is that the masses of the extra Q, Q states are likely to fall within the bounds imposed by 
string unification. 

Contrary to popular belief, string models can be tested. Moreover, if models which 
satisfy all known phenomenological constraints can ever be found, then they will become 
candidates for the theory of everything until they are defeated by future tests or until 
string theory hands down its final verdict, whichever comes first. 
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Table I: The rank and determinant [^4| of the Cartan matrices C of all the Lie algebras. 
As usual, Ai = SU(l + 1), B { = SO {21 + 1),Q = Sp(2l), A = SO {21). 



Algebra 


detC 


rank 


M 


l + l 


I 


Bi 


2 


I 


Q 


2 


I 


Di 


4 


I 


Eq 


3 


6 


Ei 


2 


7 


Eg 


1 


8 


F 4 


1 


4 


G 2 


1 


2 



Table II: The set of spin-structure basis vectors which are common in the construction 
of all our models. The first entry corresponds to the complexified ijj^ and the next 18 to 
the six left-moving triplets (x^, y^, a/). The first 12 right-moving entries (to the right of 
the colon) correspond to the real fermions y , u , and the last 16 correspond to complex 
fermions. A 1 (0) stands for periodic (antiperiodic) boundary conditions. We also use the 
symbols 1 8 = 11111111 and 8 = 00000000. 



1 = 


(1 111 111 111 111 111 111 


111111 111111 11111 111 i 8 ) 


s = 


(1 100 100 100 100 100 100 


000000 000000 00000 000 8 ) 


h = 


(1 100 100 010 010 010 010 


001111 000000 11111 100 o 8 ) 


b 2 = 


(1 010 010 100 100 001 001 


110000 000011 11111 010 o 8 ) 


bs = 


(1 001 001 001 001 100 100 


000000 111100 11111 001 o 8 ) 


6 4 = 


(1 100 100 010 001 001 010 


001001 000110 11111 100 o 8 ) 


h = 


(1 001 010 100 100 001 010 


010001 100010 11111 010 o 8 ) 
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Table III: Classification of a-vectors according to the number and type of residual U a ,b,c 
gauge symmetries which they break or do not break. U a : w 2 w 3 , Ub '■ y 1 ^ 6 , and U c : y 4 y 5 . 



Class 


rank 


Ua 


Ub 




1 


19 




V 


V 


2a 


18 


V 


V 


X 


2b 


18 


V 


X 


V 


2c 


18 


X 


V 


V 


3a 


17 


V 


X 


X 


3b 


17 


X 


V 


X 


3c 


17 


X 


X 


V 


4 


16 


X 


X 


X 
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Table IV: The massless spectrum of the selected 5/2 model with gauge group SU(5) x 
U(l)~ x SO(10)h x SU{A)h x U(l) 5 and 2/3/2 Yukawa set. The transformation properties 
of the observable sector fields under SU(5) x U(l)~ are as follows: F (10, 1/2), / (5, —3/2), 
l c (1, 5/2), h (5, —1). The hidden sector fields transform under SO(10)h><SU (4)^ as follows: 
T (10, 1), D (1, 6), F (1, 4). The F u Fj fields carry ±1/2 electric charges. 



Observable Sector : F (-±, 0, 0, -\, 0) Fi (-±, 0, 0, §, 0) 

F 2 (0,-^,0,0,0) h (0,-|,0,0,0) / 2 c (0,-i, 0,0,0) 

F 3 (0,0,±,0,-±) / 3 (0,0,i,0,i) (0,0,|,0, |) 

F 4 (-|,0, 0,0,0) F 4 (±,0,0,0,0) 

F 5 (0,|,0,0,0) / 5 (0,-|,0, 0,0) / 5 C (0,-|, 0,0,0) 

hi (1,0,0,0,0) h x (-1,0,0,0,0) 

h 2 (0,1, 0,0,0) h 2 (0,-1, 0,0,0) 

/i 3 (0,0, 1,0,0) h 3 (0,0, -1,0,0) 

h45 0,0,0) ^45 (i |, 0,0,0) 

Singlets: $ i2 (-1, 1, 0, 0, 0) $ i2 (1, -1, 0, 0, 0) 
$ 23 (0,-l,l,0,0) $23(0,1,-1,0,0) 
$3i(l,0,-l,0,0) $3i (-1,0, 1,0,0) 

045(|,|,1,O,O) 045 H, -i -1,0,0) 

0+(|,-±, 0,0,1) 0+(-±,§, 0,0,-1) 
0- (§,-§,0,0,-1) 0-(-|, 1,0,0,1) 

03,4 (J, -|,0, 0,0) 03,4 (-§,§, 0,0,0) 

7/i, 2 (0,0,0,l,0) 771,2 (0,0, 0,-1,0) $0,1,3,5(0,0,0,0,0) 

Hidden Sector: 

Ti(-i,0,i,0,0) T 2 (-±,-±,0,0,-§) T 3 (-|,0,i,0,0) 
Z?i (0, -i, i, i,0) £) 2 (0, -§,1,-1,0) D 3 (-|,0,|,0,0) 
D 4 (-i-|,0,0,|) £> 5 (0,-!,|,0,0) L> 6 (0,1, -1,0,0) L> 7 (§, 0, 0, 0) 
77+f_I I _I n -I^i p- I _I n l\ 77+ _I _I n L\ 

1 1 V 4, 4, 4' u , 2/ 2 V4' 4, 4' u ,2/ 3 V4' 4' 4, U , 2 J 

^4 (— |, |, I, 0, 0) -F 5 + ( — |, |, — |, 0, 1) F 6 + ( — 1, |, — 1, 0, — 1) 
fr~ (—1 III (—1 i i _i /7~ /'I _i i () 

- 17 1 \ 4,4,4,2' 2' 2 \ 4,4,4, 2' 2' 3 V4, 4 ,4, u , 2' 

^4 ( _ i> i> i,o, -^) f 5 (-1,-1,1,0,-5) f 6 (-|, |, _ i,o,o) 
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